Collective Two-Atom Effects and Trapping States in the Micromaser 
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We investigate signals of trapping states in the micromaser system in terms of the average number 
(n) of cavity photons as well as a suitably defined correlation length of atoms leaving the cavity. In 
the description of collective two-atom effects we allow the mean number, e, of pump atoms inside 
the cavity during the characteristic atomic cavity transit time to be as large as of order one. The 
master equation we consider, which describes the micromaser including collective two-atom effects, 
still exhibits trapping states for e 7^ 0, even for a mean number of atoms inside the cavity close 
to one. We, however, argue more importantly that the trapping states are more pronounced in 
terms of the correlation length as compared to (n), i.e. we suggest that trapping states can be more 
clearly revealed experimentally in terms of the atom correlation length. For axion detection in the 
micromaser this observable may therefore be an essential ingredient. 

PACS numbers 32.80.-t, 42.50.-p, 42.50.Ct 



INTRODUCTION 



The micromaser 0, 13 belongs to one of the simplest 
systems in cavity quantum electrodynamics. For this rea- 
son plenty of experimental as well as theoretical investi- 
gations have been carried out. A micromaser consists of 
a beam of identical two-level atoms passing through a 
cavity and interacting with a single mode of the electro- 
magnetic field. Usually the atomic flux is so small that 
at most one atom interacts with the cavity mode simul- 
taneously. This is the situation in the standard one-atom 
micromaser theory. For larger values of the atomic flux, 
and also due to statistical distribution of atomic velocities 
of the pump atoms, in an actual experimental situation 
there is, however, a finite probability for two or more 
atoms to be simultaneously present in the cavity. An 
extension of the standard micromaser theory to include 
collective effects of many atoms is therefore needed. This 
has e.g. been done in Refs.3ll|l|l|l|l[l3, where two- 
atom effects have been taken into account. Investigations 
of collective effects have been simulated by their signifi- 
cance on the mean photon number in the cavity field. At 
zero temperature, so called trapping states of the cav- 
ity radiation field 0, [T]| may appear in the micromaser 
system . These trapping states corresponds to a 

truncation of the steady-state photon distribution, and 
thereby to the appearance of a minimum in the mean 
photon number as a function of the atom cavity transit 
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time. The preparation of such trapping states, in particu- 
lar the vacuum trapping state, is essential in the design of 
axion detectors using the micromaser system (see Ref. [T^ 
and references cited therein). It has been argued that the 
appearance of trapping states maybe suppressed when in- 
cluding collective atom effects Recently, there 
has, however, been some disagreement on the degree of 
this suppression 

Several effects in the micromaser can affect the system 
dynamics and destabilize the trapping states by increas- 
ing the number of the cavity photons. The trapping effect 
in terms of the mean photon number is therefore difficult 
to observe experimentally. In the present paper we will 
therefore not only consider the photon number but also 
a suitably defined correlation length of the atoms leaving 
the micromaser. 

The paper is organized as follows. In Section|n]we out- 
line the dynamics of the standard one-atom micromaser, 
i.e. we present the well known master equation for the 
photon field in this case. In the following Section we ex- 
tend this theory to the two-atom situation, with a finite 
probability of having two atoms simultaneously in the 
cavity. The numerical results are presented in Section llVl 
together with a discussion. 



II. ONE-ATOM MICROMASER DYNAMICS 



The continuous time formulation of the dynamics of 
the one-atom micromaser is very well known (see e.g. 
Refs-OIH). We now outline a typical realization of the 
one-atom micromaser system. The pump atoms which 
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enter the cavity are at resonance with the radiation field 
of the cavity. The vector p formed by the diagonal density 
matrix elements of the photon field then obeys the master 
equation 



dp 
'dt 



-7£ip 



where 



(1) 



(2) 



Here the dimensionless atomic fiux parameter is = 
-R/7, where R is the rate of injected atoms and 7 is the 
damping rate of the cavity. The cavity damping is de- 
scribed by the matrix 



{Lc)nm = {nb + l)[n5„ 
+ nb[{n+ \)5„ 



(n + l)6n+l^ra] (3) 



where the one- atom operator matrix lAi is given by 



m-t-1 



(4) 



provided that all the pump atoms are prepared in the 
excited state. Here we have defined qn ~ s\VL'{gT \/n) — 
siv? ( 9 yJn/N ) and introduced the natural dimensionless 



pump parameter 6 = gT\J N in terms of the atomic tran- 
sit time T. Furthermore, g is the single photon Rabi fre- 
quency at zero detuning of the Jaynes-Cummings-model 
[l7|. Due to conservation of probability, the matrices 
above satisfy 



{Lc)nm — 



n=0 



and 



(5) 



(6) 



n=0 



for any given integer m > 0, as can be verified analyti- 
cally. The stationary solution of Eq. Q is then given by 
II ll 



Po 



nuy, 
. (1 



rib m + A^(7„ 



(1 + rib) m 



(7) 



The overall constant po is determined by X^^o^*" ~ 
The mean number of photons in the cavity is given by 



III. TWO-ATOM MICROMASER DYNAMICS 



Since the average number of pump atoms inside the 
cavity is e = ri? = "f^fN9/g^ this parameter naturally 
parameterize the probability of collective pump atom ef- 
fects. In the present paper, we are not restricted to e <C 1 
as in Refs. 0, , and we will even consider a case with 
e w 1. 

Let us assume that the atoms arrive at the micro- 
maser cavity with Poissonian arrival statistics. The prob- 
ability that the distance between two atoms exceeds the 
length of the micromaser cavity is then exp (— e) (see e.g. 
Ref. 4]). A one-atom event occurs whenever an atom 
is separated by at least the cavity length from both its 
neighbours in the beam. Likewise, a two-atom event oc- 
curs whenever two atoms are separated from each other 
by less than the cavity length and the nearest neighbours 
of the pair are further away then the cavity length. The 
probability that a randomly-chosen atom in the beam is 
a n-atom event then is l4| 



Pn=e- 



(8) 



The probabilities in Eq. (|SJ) are normalized in such a way 
that the total atomic fiux R is fixed, i.e. X^J^i nRPn = 
R. A generalization of the one-atom master equation to 
one that includes one-atom events as well as two-atom 
events has then been proposed in Ref.Q, i.e. 



dp 
'di 



where 



£■2 = Lc - N Pi (Ui - 1) - N P2 {U2 - l \ 



(9) 



(10) 



Here the modified one- and two-atom probabilities. Pi 
and P2, are given by ^] 



Pi - 
P2 = 



Pi 



Pi + 2P2 

P2 
Pi + 2P2 



(11) 
(12) 



normalized in such a way that Pi + 2P2 = 1. The one- 
atom operator Ui is given by Eq.Q). In particular, for 
Pi — 1 and P2 = the two-atom generator £2 is reduced 
to the one-atom generator as defined in Eq. Q . The two- 
atom operator U2 is defined in such a way that U2 Pf (t) 
gives the reduced density operator for the radiation field 
for a two-atom event H, i.e. 



1^2 P fit) 



ds w{s) 

X Trat,at2{"i(al2,s)M2(T 
X Pati ® Pats ® /'/W } ■ 



(13) 



I ui(ati, s) 
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Here Ui(ati,s), for i = 1,2, denotes the combined one- 
atom (ati) and radiation field unitary time-evolution op- 
erator for a time- interval s. Similarly, U2{t — s) denotes 
the combined two-atom- and radiation field unitary time- 
evolution operator for the time- interval t — s. Tracing 
over the atomic states in Ea. H13(l gives the reduced den- 
sity operator of the radiation field for a two- atom event. 
The integrand in Eg. ljl^ll describes a two-atom event, in 
which case the two atoms arrive a time s apart and over- 
lap for a time t — s during which both atoms are present 
in the cavity. Since we assume that the arrivals of the 
atoms in the cavity are uncorrelated events of a Poisso- 
nian process, the two-atom events enters in an averaged 
form with the conditional probability density distribution 



w{s) = R 



(14) 



i.e. the Poisson probability for two atoms arriving a time 
s apart conditioned by an overlap during the atomic tran- 
sit time T. The distribution w{s) is therefore normalized 
to one in the atomic transit time r. 

In the quantum trajectory approach to collective two- 
atom effects 6, 7], an ensemble average is calculated, 
where the overlap s is randomly selected according to the 
probability distribution Ea. lll4l) . Provided the system is 
ergodic, such an ensemble average should corre spo nd to a 
time-average according to Ea. (|13|l . In Refs. 5, 10] a sim- 
ilar time-average has been considered. Recently it has 
been debated to what extent the micro-maser system is 
ergodic when multi-atom effects are taken into account 
[a ISj , at least when considering the average number of 
cavity photons number (n) . Trapping states appear to be 
more suppressed in the quantum trajectory approach 
as compared to the time-average approach j3| as defined 
by Ea. (|13|l . Below we will, however, consider an alter- 
native observable, which is more clearly sensitive to the 
presence of trapping states than (n) but which, at the 
same time, is not very sensitive to small changes in the 
physical parameters of the micromaser system. For this 
reason we conjecture that calculations of this observable 
in the quantum trajectory approach and the time-average 
approach will lead to similar results. In the present paper 
we will restrict ourselves to the conventional time-average 
approach as in Ea. (|13|l . 

Explicit and correct expressions for the evolution op- 
erators ui and U2 can be found using Ref . |^ . A straight- 
forward calculation then gives 



{U2)nm. ^ ds w{s) (Z^2(s))„m , (15) 

Jo 

where we make use of the notation 



ms))r. 



(16) 



In the expression for the matrix element (U2{s))nm we 
have defined 



and 



/ ~ ~ 1 2 

2 v 2€7T, C-nSn ^n+l \ ; 

Kab,n = I Sn+1 C„+l [ ( 1 + 2e„ ) - 1] 

+ '\/2en c„s„ [ s„+i — c„_|_]^ ] I 

+ \-[cisn+l + V2e^Cn-lSnCn+l]Cn+2 



(17) 



+ [ 2^e„(l - e„) si Cn+i 

, (18) 



as well as 



e?i) C,i_lS„+l Sn+l ] Cn+2 + 

I • (19) 



In these expression we make use of the convenient defi- 
nitions 



1 



2n + 3 



(20) 



and trigonometric functions for the one-atom events 

Sn = sin[gsVri] , (21) 
c„ = cos[gs\/n] , (22) 

as well as for the two-atom events, i.e. 



sin[5(T — s) \/n-\- 3/2] 
cos[.g(r - s) + 3/2] 



(23) 
(24) 



The expressions we obtain are the same as in Ref. Q • We 
also mention that, due to unitarity, the two- atom matrix 
satisfies 



(^2)nm — 1 



(25) 



n=0 



for any given integer m > 0. This can be verified an- 
alytically. It is also now straightforward to verify the 
relation U\ = U2 (t) , which corresponds to the case when 
the atoms in the beam are separated in time by s = r, 
i.e. no overlap. 

In particular, let us consider the limit in which the 
mean number of atoms in the cavity is small, i.e. e <C 1. 
In this case the micromaser is close to the one- atom 
maser situation. From Eqs. Hll(l and H12|l we then ob- 
tain Pi « 1 — 2e and P2 ~ e. Furthermore, the condi- 
tional Poissonian probability distribution Eq. I|14|l is then 
reduced to w^s) ~ 1/t. We then obtain 
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dp 
'dt 



« -7 ic-iV(l-2e)(Z^i-l) (26) 0.8 



Ne{U2-l)\p, 



where 



ds {U2{s))r 



(27) 



The master equation E g J26(l appears to be of the same 
form as the result in Ref.0(see also R.efs. flfllT^ I. There 
is, however, one important difference. In Ref.P| it is 
assumed that both e and the one-atom generator lAi — \ 
are to be treated as small entities. This e.g. means that 
U2{t) = l + 2{Ui-\) + 0{{Ui-lf). In order to compare 
Ea. (PS|l with the corresponding expression of Ref.0, we 
should therefore in general only retain terms to at most 
second order in the one-atom trigonometric functions c„ 
and s„ in Eas. lfT7|) - ((T^ . With these approximations we 
arrive at the same expression as in Ref. Q as is also used 
in Ref. ^3 ■ our numerical simulations to be presented 
below we make use of the exact expression for . 

Our final expression for £2 as given in Ea. H10|l is the 
same as in Ref. . It is rather straightforward to imple- 
ment the matrix elements of C2 in a numerical routine 
in order to find the stationary probability eigenvector p 
corresponding to the eigenvalue Ao = of £2 ■ This eigen- 
vector p is to be used in evaluating the mean number of 
photons (n) in the cavity. The next-to-lowest eigenvalue 
Ai of £27 which we find numerically, will then determine 
typical time scale for the approach to the stationary sit- 
uation as described by the eigenvector p. The joint prob- 
ability for observing two atoms, with a time-delay t be- 
tween them, can now be used in order to define a corre- 
lation length ^^{t) as first suggested by Lautrup et al. 
[3 (for a recent discussion see Ref. ilOj ) . At large times 
t 00, the atomic beam correlation length is then given 
by by El, i.e. 



lA{t) 



-t/iA 



(28) 



which then is determined by Ai, i.e. = 1/Ai. In 

passing we notice that for photons a similar correlation 
length was also defined in Ref. '16] in terms of cor- 
relations of the number of photons. It was shown that, 
quit naturally, these correlation lengths are identical, i.e. 

= = C 0- tbe numerical work presented be- 
low, it turns out to be sufficient to use 200x200 matrices 
for £2 in order to obtain the accuracy of the graphs as 
presented in the present paper. 

IV. RESULTS 
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FIG. 1: The order parameter (x) = (n/N) as a function of 
for TV — 50, ni, = 0. The upper graph shows the order 
parameter for the one-atom micromaser, i.e. e = 0. The 
lower graph shows the micromaser when two-atom events are 
included with e = 0.05, 0.1, 0.4. 



It was argued in Refs.[^0| and in Ref. 7] that trapping 
effects, which occur at 9 — kn^J N/n, {k,n = 1,2,...), 
may be significantly suppressed due to the inclusion of 
collective two-atom events. For reasons of comparison, 
FigQand Fig|21in the present paper are plotted for the 
same parameters as in Refs.p, [l^- In Fig ^ we show 
the results of a numerical evaluation of the order param- 
eter (2:), where x = n/N and ( ) as above denotes an 
average with respect to the stationary solution p when 
two-atom collective effects are included. For e = 0, the 
non-trivial and detailed structure of the various trapping 
signals , i.e . the observed dips in (x), has been explained 
in Ref.|l9j|. The numerical precision of our calculations, 
as well as in Ref.^l, is sufficient to reveal the exact 
large- TV structure of these dips in the order parameter 
{x). This is not the case in Ref.Q. We, furthermore, 
observe a discrepancy between our numerical results and 
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FIG. 2: The correlation length 7^ as a function of 9 for — 
50, n;, = 0. The upper graph shows the correlation length for 
the one-atom micromaser, i.e. e = 0. The lower graph shows 
the micromaser when two-atom events are included with e — 
0.05,0.1,0.4. 



the findings of Ref . |3] . 

As mentioned above, it may, however, be difficult to 
experimentally discriminate between the collective effects 
and other corrections, like the effects of detection efficien- 
cies, when only considering the order parameter (x). As 
seen from Fig|21 where the correlation length 7^is plot- 
ted on a logarithmic scale, trapping states are more pro- 
nounced in terms of this observable as compared to the 
order parameter (x) as given in Fig^ Even for such a low 
value of the atomic flux parameter as iV = 10, and with 
rib = 0.001 as in FigOl the trapping states are clearly 
revealed. These results therefor convincingly shows that 
trapping states can be even more clearly revealed exper- 
imentally in terms of the correlation length 7^ as com- 
pared to the order parameter (x). In this context it is 
important to notice that the effect of detection efficien- 
cies on the correlation length 7^ easily can be taken into 
account in terms of a simple renormalization of the ex- 
perimental data as discussed in detail in Ref. ,10l| . 
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FIG. 3: The correlation length 7^ as a function of 6 for 
various values of e. The vertical lines indicate the trappings 
values at 61 = 7r//V(l/V3, 1/^, 1, 2/^3, ^2, ^,2). 



the corresponding results in Refs.[^0|. This discrep- 
ancy is mainly due to fact that we make use of the exact 
form of the two-atom generator which was not the 
case m Refs.ilia. In fact, it turns out that the exact 
expression for the conditional probability density Eq. I|14|) 
and the modified probabilities Ea. (|lll) and Ea. l|12() are 
of less importance. We have verified this by reconsid- 
ering the numerical evaluations of Ref. in terms of 
the Eas. ((n|l . ^fTT}i and (fT^ instead of using their small e- 
expansions, and we find only small numerical corrections 
to results presented in Ref. 10] . An exact treatment of 
the two-atom generator U2 therefore, remarkably, leads 
to smaller collective effects as compared to the results of 
Refs.0, 01 . Most of the dips in (x) corresponding to 
trapping states actually do not vanish for e 7^ in our 
case, even for mean number of atoms inside the cavity 
as large as e = 0.4. This conclusion is in agreement with 



In the experimental study of trapping states in Ref. |l2l 
the parameters are varied in such a way that the 
parameter e is not constant but is given by e = Rt = 
In Fig^ we have chosen the cavity temper- 
ature to be T = 0.3 K, corresponding to rib ~ 0.054, 
and a low pump rate corresponding to = 10. These 
particular choices of parameters are the same as in Wei- 
dinger et al. |12| |. in which case two-atom events are rare, 
i.e. e < e{0 = 20) = 0.016 < 1. The solid line in 
FigElcorresponds to e = 0, i.e. no collective effects taken 
into account. Comparison with the dashed line, where 
two- atom collective effects are taken into account, gives 
basically the same result. This is not unreasonable as the 
mean number of photons inside the cavity is very small 
in this case. We also mention that, for the particular pa- 
rameters under consideration, the result in Ref. [13] (see 
their Fig. 6) gives basically the same result. 
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FIG. 4: The order parameter (a;) and the correlation length 
7^/5 as a function of 6 with no collective effects taken into 
account (e = 0) or e = Rt = O^/N'y/g. The Rabi fre- 
quency is as in the experimental study of trapping-states in 
Ref-Eal, i.e. g = 39kHz and 7 = 10 s~\ The vertical 
hnes indicate the experimentally observed trappings values at 
e = nVN{l/V3,l/V2, 2/V6, 2/^5, 1, 2/VS, v^) as reported 
in Ref.fil. 
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FIG. 5: The correlation length 7^ as a function of 6 for 
various values of e. The vertical lines indicate the 9 values of 
the large A'^ one-atom maser phase transitions at So = 1, 6^01 ~ 
6.6610,6*^2 ~ 12.035,013 ~ 17.413. With (x) as an order 
parameter, the transition at 9 = 9n is second-order while the 
other are first-order transitions figll. 



find that collective two-atom effects are small. 



In R,efs.[T^ IT^ the phase structure of the micromaser 
system in terms of the physical parameters at hand is in- 
vestigated. It is shown that there are two types of transi- 
tions, a thermal-to-maser transition and maser-to-maser 
transitions. In the large N limit, the thermal-to-maser 
transition occurs at 9q and the maser-to-maser transi- 
tions occur at d^i^^i for fc > in the absence of collective 
effects. For given values of the parameters, these values 
of the pump parameter are computed numerically. At 
finite N and including collective pump atom effects, we 
can see from e.g. Fig|2and in particular Fig[Slthat sig- 
nals of the large N phase transitions in the correlation 
length 7^ are still clearly exhibited close to the corre- 
sponding one-atom maser phase transitions. As pointed 
out in Refs. 0, 0, 0| , the critical point 9q of the first 
second-order maser transition is clearly exhibited in both 
(x) as well as log(7^). This fact is illustrated for 7^ in 
Fig|Sl Furthermore, in Ref.[To| using a perturbative rea- 
soning, it is argued that, for sufficiently small e, the crit- 
ical parameters B^^i oc y/N of first-order transitions are 
changed to O^^.^-^ = exp(-e)0|;j,^i and that log(7^)crjt 
is changed to log(70crit = exp(-e) log(7^(6' = Ol^+i))- 
This change of the correlation length was also in agree- 
ment with the numerical result of Ref . |lOl | . When using 
the exact expression for the two-atom generator these 
perturbative considerations are, however, not valid as e.g. 
can be seen from our numerical results presented in Fig|Sl 
In this figure we see that the peak values of log(7^)crit are 
actually shifted slightly to higher rather then to smaller 
values of the pump parameter 0. In Fig[Sl we have also 
evaluated the correlation length for e = 1 and we still 



V. CONCLUSIONS 



In conclusion, we have studied two-atom collective ef- 
fects in the micromaser system. We have confirmed the 
result of Ref. S| that the characteristic effects of trapping 
states on the order parameter are not as suppressed as 
was previously suggested 0, 0| . We have seen that the 
major reason for this is connected to the exact treatment 
of the two- atom generator. More importantly, the pres- 
ence of trapping states are more pronounced in terms of 
the correlation length 7^ as compared the order param- 
eter (x). This conclusion remains true when collective 
two-atom effects are taken into account even for a large 
average value of atoms inside the cavity. As seen in Fig01 
which refers to an actual experimental situation 12|, the 
vacuum trapping state, which for the actual micromaser 
parameters occurs at 9 = tt^/N !v 9.94, is clearly more 
visible in the correlation length than in the average num- 
ber of photons. This conclusion also remains true when 
collective two-atom effects are taken into account. An 
experimental preparation of the vacuum trapping state 
is essential when considering axion detection in the mi- 
cromaser system JJj. We conclude that trapping states 
can be more clearly revealed experimentally in terms of 
the atomic correlation length than in terms of the natu- 
ral micromaser order parameter (n/N). This observation 
is a ppa rently in accordance with the recent analysis of 
Ref. Ill . 
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